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ABSTRACT 

The stationary spherically symmetric accretion flow in the Schwarzschild metric has 
been set up as an autonomous first-order dynamical system, and it has been studied 
completely analytically. Of the three possible critical points in the flow, the one that 
is physically realistic behaves like the saddle point of the standard Bondi accretion 
problem. One of the two remaining critical points exhibits the strange mathematical 
behaviour of being either a saddle point or a centre- type point, depending on the values 
of the flow parameters. The third critical point is always unphysical and behaves like 
a centre-type point. The treatment has been extended to pseudo-Schwarzschild flows 
for comparison with the general relativistic analysis. 
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1 INTRODUCTION 

To researchers in astrophysics and general relativity, phys- 
ical models of spherical symmetry have an abiding appeal. 
One especial advantage with these models is that almost 
always they lend themselves to an exact mathematical anal- 
ysis, and in the process they allow a very clear insight to be 
had into the underlying physical principles. For this reason 
in particular, spherically symmetric models frequently serve 
as a firm foundation for constructing theoretical models of 
physical systems involving more realistic and, therefore, un- 
avoidably complicated features. 

Studies in accretion are no exception t o this practice . 
Ever since the seminal work published by iBondil (|l952t ). 
that effectively launched the subject in the form in which 
it is recognised today, the problem of spherically symmetric 
flows has be en revisited time and agai n from various 
angles |Parkeri, 1958,., 1966: Axford fc Newman,, 1967.: .Balazs 
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The original work of iBondil (|l952t ) introduced formal 
fluid dynamical equations in the Newtonian construct of 
space and time to study the stationary accretion problem. 
From here the connection to the proper general relativis- 
tic framework was n ot too long in coming. In particular 
it was iMichell l| 197^ 1 who made an important early foray 
into the general relativistic domain, which was followed by 
a spate of other works, o f which some (to mention a few), 
along with the paper of iMichell (Il972l ). addressed various 
aspects related to the critical b ehaviour of general relativis- 
tic flows in spherical symmetry ("Begel manll 19781 : IBrinkmannI 
|l980; Malcc 1999; Das fc Sarkar 2001). 

The work being presented here is also along the same 
lines. Its direct purpose is to construct a pedagogical theory 
to understand the nature of the critical points of station- 
ary spherically symmetric flows in the Schwarzschild met- 
ric, after starting with the set of basic stationary equations 
which govern the flow. To the extent that critical solutions 
— speciflcally transonic solutions in regard to spherically 
symmetric flows — can only pass through some critical 
points (which must arguably be saddle points), this treat- 
ment will also have a bearing on a very important gen- 
eral issue in accretion studies — the manner in which a 
compressible astrophysical fluid passes (either continuously 
or discontinuously) from infinity to the event horizon of 
a black hole. In addressing this question the mathemati- 
cal method that has been adopted is a dynamical systems 
analysis, which is always an effective tool for researchers in 
non-linear dynamics, and which has been tried successfully 
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before in a related astrophysical system, namely multitran- 
sonic flows i n an inviscid and thin p seudo-Schwarzschild ac- 
cretion disc (jChaudhurv et al]l2006l ). While it is to be natu- 
rally expected that the same treatment could be carried over 
directly to a rotational flow in the Schwarzschild metric, it 
has also been a worthwhile exercise to consider spherically 
symmetric flows flrst, as a suitably uncomplicated starting 
point into more intricate general relativistic problems. An 
immediate advantage in studying this relatively simple sys- 
tem has been that the mathematical treatment could be 
carried out with full analytic rigour, something, as it has 
been stressed right at the beginning, that has inspired the 
choice of the spherically symmetric model in the flrst place. 
It has been found here that under practical restrictions on 
the conditions for critical behaviour, it will be possible to 
gain a complete mathematical understanding (in the sense of 
producing final results which are absolutely non-numerical) 
of the nature of the critical points (not all of whom may be 
physically realistic), and the pattern of the solution topolo- 
gies in the neighbourhood of those points. In no way do 
these results contradict any of the earlier findings (|Miche]| 
Il972 : Bcgclman 1978), and if anything, many surprisingly 
new features in the flow, hitherto unrecognised, have been 
revealed. 

The dynamical systems treatment has also been car- 
ried out on pseudo-Schwarzschild spherical flows driven 
by so me of the establ i shed pseudo-Newtonian poten- 
tials (iPaczvnski fc Wiital Il980l : iNowak fc Wagoned Il99ll : 
lArtemova et al.l 19961 ) to check for the consistency of this 
approach with the fully general relativistic methods. All the 
results have been in complete qualitative compatibility with 
one another in the sense that there is only one physically 
meaningful critical point in the flow through which a so- 
lution could pass transonically, connecting inflnity to the 
event horizon of the black hole. In one important detail, 
however, a quantitative difference has appeared. For accre- 
tion governed by cold ambient conditions, transonicity has 
been shown to be very much more pronounced in a properly 
relativistic flow, than in a pseudo-Schwarzschild flow. 

The case for a dynamical systems approach in study- 
ing a proper general relativistic problem has been argued 
cogently in this work. The mathematical methods demon- 
strated here can very well be applied to the more involved 
accretion disc system, described both by the Schwarzschild 
metric and the Kerr metric. This particular study will be 
reported separately. Meanwhile, following the treatment on 
pseud o-Schwarzschild flows reported by IChaudhurv et al] 
l|2006l ). the present work on spherically symmetric flows may 
be considered to be the second in a series that pedagogically 
underlines the conspicuous advantages of applying dynam- 
ical systems methods in standard astrophysical fluid flow 
problems. 



2 THE GENERAL RELATIVISTIC FLOW AND 
ITS FIXED POINTS 

In this general relativistic treatment of a spherically sym- 
metric, stationary, compressible fluid flow, the two relevant 
flow variables will be the radial inflow three-velocity, v, and 
the local proper mass density of the fluid, p. The radial co- 
ordinate of the flow is scaled by the Schwarzschild radius. 



Tg — 2GMbh/c^ (here Mbh is the mass of the black hole), 
with any characteristic velocity in the flow being scaled 
by c. Setting G = c = 1, the general relativi stic a nalo g ue 
of Bernoulli's e quation will be given as (jChakrabartiill990l . 
ll996l : [Da3[200i ') 



£ = hvt = 



p + e 



(1) 



with the pressure, p, connected to the density, p, through 
an equation of state, p = kp~' , and the speciflc enthalpy, h, 
expressed as 
p + e 



h = 



(2) 



in which the energy density, e (which includes the rest mass 
density and the internal energy), is to be set down as, 

P 



7-1 

It is then possible to arrive at the relation 



£ 



7 - 



7-1 



1 - ?;2 



(3) 



(4) 



with the speed of sound, Ca, deflned under conditions of con- 
stant entropy, 5, as 



2 dp 
oe 



(5) 



Through the equation state, p = fcp^, the speed of sound 
can be connected to the density, p, as 



(6) 



7^ (1 — nci) 



with w being given by the u sual definition of the polytropic 
index (jChandrasekharll 19391 ) as n = (7 — 1)"^. 

The stationary continuity condition, on the other hand, 
will give another relation connecting the velocity and density 
fields as. 



1 



(7) 



with m being an integration constant, which is to be physi- 
cally identified as the mass accretion rate. 

It is now easy to see that equations ([4]), and ((Tjl 
will give a complete description of the flow system. Making 
use of equation ((6} in equation ([7|, and then going back to 
equation (j4|, it will be possible to express the gradient of 
solutions in the r — plane as 



d , 2^ 



[cl (4r - 3) - 1] 



(8) 



r (r — 1) (u^ — Ca ) ' 

with Ca being used as a characteristic local scale of velocity 
in the fluid, against which the local bulk velocity of the flow, 
V, is to be measured. 

From the foregoing expression it is evident that there 
will be non-trivial singularities under the conditions r = 1 
and = Ca, unless the numerator in equation ^ vanishes 
simultaneously. The condition r = 1, of course, corresponds 
to the behaviour of the flow on the actual event horizon of 
the black hole (where = 1), but of immediate interest 
is the condition — c^, which will give a critical condi- 
tion for the flow at r > 1, if and only if the requirement of 
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Cs (4f — 3) = 1, is simultaneously satisfied in the numerator. 
In that event the critical point conditions in the flow will be 
expressed as 



1 



4rc 



(9) 



with the subscript "c" labelling the critical point values. 

The next logical step from here is to represent the crit- 
ical point coordinates in terms of the parameters of the sys- 
tem. This can be done in two ways — either by substituting 
Wc and Csc in equation Q and expressing rc as a function of 
£ and 7, or by substituting the same critical point values in 
equation ((T)) and expressing as a function of m and 7. Ei- 
ther approach is entirely equivalent to the other, and here for 
simplicity of algebraic manipulations, the former approach 
is being adopted. This will deliver a cubic equation in rc as 



0, 



(10) 



where 



^0 = 



27 



64 (£:2-l)' 



A^ = 



and 



A2 = 



(2-37) -27(7- 1) 
16 (£2 - 1) (7-1)^ 



2 (2 - 37) + 9 (7 - 1) 



4(22 -1) (7- 1) • 

The solutions of equation (|10|l can be found completely 
analytically by employing the Cardano-Tartaglia-del Ferro 
method for solving cubic equations. To that end it should 
be first convenient to define 



El 
S2 



3^1 - Al 



9A1A2 - 27 Ao - 2Al 
54 



and 



6 = (S2 



1/3 



following which, the three roots of rc can ultimately be set 
down as 



rci 



' 3 
A2 



+ (6+6), 



r-c2 = -— -(6+6) + i^ (6-6) 



and 

rc3 = 



^^-(6+6)-i^ (6-6): 



(11) 

(12) 
(13) 



respectively. 

The sign of ^ should be crucial in determining the na- 
ture of the roots. If > 0, then only one root of rc will be 
real, while for < 0, there will be three real roots. Since it 
is obvious that will have a dependence on E and 7, it will 
be instructive to consider the appropriate ranges of values 
for these two parameters. 



The parameter £ is scaled in terms of the rest mass en- 
ergy and it includes the rest mass energy itself. So it might 
be argued that a lower limit of E would be £" = 1. On the 
other hand, although E can in principle assume any value 
greater than unity, values of f > 2 will imply extremely hot 
conditions at the outer boundary, with the thermal energy 
being much greater than the rest mass energy. Such a situ- 
ation could not conceivably prevail in realistic astrophysical 
systems, and so the practically admissible range of E will be 
restricted to 1 < f < 2. The non-relativistic range of £, on 
the other hand, will be < f < 1, without any involvement 
of the rest mass energy. This range will be considered for 
the pseudo-Schwarzschild treatment in Section |4l 

The parameter 7 is likewise restricted by the range 
1 < 7 < 2. The lower limit, i.e. 7=1, corresponds to op- 
tically thin, isothermal accretion, while values of 7 > 2 will 
involve magneto-hydrodynamics in the general relativistic 
theory and the anisotropic nature of pressure. For most re- 
alistic purposes, it should be noted, the range of 7 actually 
varies from 4/3 to 5/3. Detailed discuss ions devoted to these 
issues are to be found in the literature llBegelmarJI 19781 : iDasI 
|l)04ii). 

And so it is that with the ranges of 1 < f < 2 and 
1 < 7 < 2, it should be easy to show that ^ would always be 
negative. This will consequently imply that the three roots 
of rc, as given by equations (lllf) . (|12|l and (|13|) . are always 
real, and they can be represented in terms of a new variable, 

/ \ 

B = arccos 



A2 



icj — — — + 2\/— El cos 



e + 27r (i - 1) 



(14) 



with the label j taking the values {j — 1, 2, 3}, respectively, 
for the three distinct roots. Of these roots, rc2 is always neg- 
ative and, therefore, is not of much physical interest. The 
other two roots are always positive, and of these two, the 
one at rd is always to be found at distances greater than 
the event horizon of the black hole, i.e. rd > 1. A physically 
meaningful transonic inflow solution, connecting infinity to 
the event horizon, seems to prefer the critical point at rd to 
the one at rcs, even when rcs > 1. Through the latter point, 
the flow exhibits non-physical properties like the matter in- 
flow rate, as given by equations (|6} and (O , being reduced to 
an imaginary quantity. It was also pointed out by Das ( 200^) 
that a flow associated with this point becomes superluminal 
much before reaching the event horizon of the black hole. 
Nevertheless, for all its apparent barrenness from a physical 
perspective, this critical point is not entirely devoid of some 
very interesting mathematical properties, when a dynamical 
systems approach is made to study the nature of the critical 
points of the flow. This issue will be taken up in Section |3] 



3 PROPERTIES OF THE FIXED POINTS : AN 
AUTONOMOUS DYNAMICAL SYSTEM 

So far the flow variables have been ascertained only at the 
critical points. Since the flow equations are in general non- 
linear differential equations, short of carrying out a numer- 
ical integration, there is no completely rigorous analytical 
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prescription for solving these differential equations to deter- 
mine the global nature of the flow variables. Nevertheless, 
some analytical headway could be made after all by tak- 
ing advantage of the fact that equation (|8|, which gives a 
complete description of the r — phase portrait of the 
flow, is an autonomous flrst-order differential equation, and 
as such, could easily be recast into the mathematical form 
X = X{x, y) and ij = Y(x, y), which is that of the very famil- 
iar c oupled flrst-order dynamical system ijjordan fc Smith! 
1 19991 ). 

With the adoption of this line of attack, equation ([8]) 
may be decomposed in terms of a mathematical parameter, 
r, to read as 



d , 2. 

d7(") 
dr 
d7 



(l - v') [c! (4r ^ 3) - l] 
r (r - 1) (t;^ - Cg) 



(15) 



in both of which, it must be noted that the parameter r 
does not make an explicit appearance in the right hand side, 
something of an especial advantage that derives from work- 
ing with autonomous systems. This kind of parametriz ation 
is quite common in fluid dynamics (iBohr et al.lll993) . and 
in accre tion studies especially, this a p proach has been made 
before dRav &: Bhattacharied I2OO2I : lAfshordi fc Paczviiskil 
I2OO3I : IChaudhurv et ahlbooi ). A further point that has to 
be noted is that the function Cg in the right hand side of 
equations l|15p can be expressed entirely in terms of and 
r, with the help of equations (|6]) and ([7]). This will exactly 
satisfy the criterion of a flrst-order dynamical system. 

The critical points of the foregoing dynamical system, 
as equations ^ give them, have already been identified, and 
as equation (|10p indicates, they have also been fixed in terms 
of the physical fiow parameters. Beyond this stage, the next 
task would be to make a linearised approximation about the 
fixed point coordinates and extract a linear dynamical sys- 
tem out of equations (|15|l . This will give a direct way to 
establish the nature of the critical points (or fixed points), 
which will ultimately pave the way for an investigation into 
the global behaviour of the solutions in the phase portrait 
of the fiow. Indeed, not infrequently, if the fiow system is 
simple enough, with only an understanding of the features 
of its critical points, complete qualitative predictions can be 
made about the global solutions. In this regard the classi- 
cal spherically symmetric Bondi fiow , with its single criti- 
cal p oint, provides an object lesson (|Rav fc Bhattacharied 
|2002| ). 

Expanding about the fixed point values, a perturbation 
scheme of the kind v'^ = + Sv'^ , = c^^ + 5c^ and r — 
Vc + Sr is now to be applied on equations psp . and then a set 
of coupled autonomous linear equations is to be derived from 
them. While doing so, it will also be necessary to express Sc^ 
itself in terms of 5r and 5v^ , with the help of equations (O 
and (O, as 



5c| 



4rc — 3 / (5r 



1 



.(16) 



The resulting coupled linearised equations will then read as 



-BSv' + 



4c^ 



B 1 



Tc (re — 1) 



Sr 




Figure 1. Variation of D,^ (associated with a physical saddle 
point) with respect to the parameters £ and 7 for the fully gen- 
eral relativistic spherically symmetric flow. All values of are 
positive for the chosen ranges of £ and 7. For small values of these 
two parameters, the saddle-type feature is very robust. 



d7(^^) = 



^c(^c-l)(l + ^) 



5v + B Sr, 



(17) 



with B — [j ~ 1 — Cst.)/2. Using solutions of the type Sv^ ~ 
exp(nr) and 5r ~ exp(nr) in equations (|17p . the eigenvalues 
of the stability matrix associated with the critical points will 
be derived as 



+ (2rc -C) 



3(rc- 1) 



■ 7 



4rc 



(18) 



where 



C = 



4(7-l)r-c~(37-2) 
4rc - 3 



Once the position of a critical point, r^, has become 
known, it is then quite easy to determine the nature of that 
critical point by using rc in equation (|18p . Since is a 
function of £ and 7, it effectively implies that 57^ can, in 
principle, be regarded as a function of the fiow parameters. 
From the form of fi'^ in equation (|18|l . a generic conclusion 
that can be immediately drawn is that the only possible 
critical points will be saddle points and centre-type points, 
and for the former, > 0, while for the latter, < 0. 
This is entirely to be expected because the physical system 
under study here is a conservative system, very much like, 
by analogy, the undamped simple harmonic oscillator, the 
fixed point s of whose phase portr ait also manifest identical 
properties (| Jordan fc Smithlll999^ . 

Of the three critical points, as implied by equation (|14p . 
the one given by the unphysical negative root, rc2, is always 
a centre- type point. This is something that by itself is also 
an equally unphysical trait as far as transonic accretion is 
concerned, where the whole objective is to have a solution 
that will connect infinity to the event horizon of the black 
hole, and in doing so will cross the sonic barrier with a finite 
gradient. 

This requirement, on the other hand, is very eminently 
met at the critical point fixed at rd, which is always a sad- 
dle point and allows a physical transonic solution to pass 
throught itself without any hindrance whatsoever. The be- 
haviour of this critical point has been graphically depicted 
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Figure 2. The general relativistic flow has a strange critical point 
which behaves sometimes like a saddle {Q? > 0) and at other 
times like a centre-type point {Q? < 0), depending on the values 
of the parameters E and 7. The dependence of does not exhibit 
widespread deviations, except for small values of E and 7. 



in Fig. [T] An interesting fact that emerges from the plot 
is that for low values of £ and 7, there is a strong growth 
behaviour for . In fact, as £ and 7 approach unity, the 
value of increases by four orders of magnitude than what 
has been scaled along the vertical axis of Fig. and this is 
exactly how it should be. Apart from the sign of Q?' — al- 
ways positive and, therefore, indicative of a saddle point — 
its magnitude also conveys quantitative information about 
the strength of the saddle-type behaviour; the greater its 
magnitude, the more prominent in effect will be the tran- 
sonic behaviour of the flow (which will pass through the 
saddle point). When £ and 7 approach unity, it will cor- 
respond more closely to a cold, isothermal distribution of 
matter, and this will make the flow easily submit to the 
strong gravitational influence of the black hole (a hotter and 
more pressure-dominated flow will be capable of building a 
much greater resistance against gravity). Transonicity can 
only occur when gravity triumphs over all other effects, and 
the general features of the flow indicated by Fig. [l] is an 
emphatic endorsement of physical transonicity in a general 
relativistic scenario. 

On the other hand, a most intriguing and counter- 
intuitive behaviour is to be encountered at the critical point 
characterised by r^^, which is placed between rd (always 
a saddle point) and rc2 (always a centre-type point). De- 
pending on the values of the parameters £ and 7, this point 
exhibits the properties of both a saddle point (fi^ > 0) and a 
centre-type point {Q? < 0). As far as the former case is con- 
cerned, this is a very curious state of affairs indeed, because 
conventional wisdom about well-behaved fixed points will 
have it that no two adjacent fi xed points can both be saddle 
points l| Jordan fc Smithll 19991 ). quite contrary to what is be- 
ing seen here — a saddle point (obviously realistic and phys- 
ically meaningful) and a centre-type point (however physi- 
cally unrealistic) fianking a point which, on some occasions 
at least, behaves like a saddle (and on other occasions like a 



^ This reduced scaling of the vertical axis will give a much better 
resolution of the overall behaviour of against the parameters 
on which it is dependent. 



centre). It is difficult to provide an analogy for this kind of 
behaviour from any other area in physics. One might con- 
jecture that this waywardness could be intimately connected 
to the divergent behaviour (like superluminal motion) exhib- 
ited by solutions associated with this critical point (when it 
behaves like a saddle point). A quantitative graphical un- 
derstanding of the nature of this critical point has been con- 
veyed in Fig. [2I Once again it is to be seen that as £ and 
7 both assume values closer to unity, strong evidence of a 
saddle-type behaviour results. 



4 THE PSEUDO-SCHWARZSCHILD 

APPROACH : A COMPARATIVE STUDY 

Frequently in studies of black hole accretion, it becomes con- 
venient to dispense completely with the rigour of general 
relativity, and instead make use of an "effective" pseudo- 
Newtonian potential that will imitate general relativistic ef- 
fects in the Newtonian construct of space and time. In that 
event the relevant stationary equations for the compressible 
spherically symmetric fiow will look like 

,J: + i^ + ^'(,)=0 (19) 
ar p ar 

and 



(20) 



Ar 

respectively, with the former being the familiar Euler's equa- 
tion and the latter the equation of continuity (jChakrabartil 
!l990l : lFrank et al.ll20o3 ). In equation (US]), <i>{r) is the gen- 
eralised pseudo-Newtonian potential driving the flow (with 
the prime denoting its spatial derivative), and p is the pres- 
sure of the flowing gas, which is related to the density by 
the usual polytropic prescription. The local speed of sound, 
with which the bulk flow will have to be scaled, is defined 
by Cg — dp /dp, following which the connection between the 
p and Cs could be established as 



7A; 



(21) 



whose form may, for academic interest, be compared with 
the general relativistic analogue given in equation (|6]). 

Making use of equation (|21|l in equation (|20|l , and then 
going back to equation (|19p . will lead to a relation for the 
gradient of solutions, which will read as 



Ar 



2v'^ [2cl - r(j)'{r)\ 
r (u^ — c|) 



(22) 



The critical points in the flow will be derived from the 
standard requirement that the flow solutions will have a fi- 
nite gradient when they will cross the sonic horizon, which 
will mean that both the numerator and the denonimator 
will have to vanish simultaneously and non-trivially. This 
can only happen when 



c0'(rc) 



(23) 



which gives the critical point conditions, with the subscript 
"c" labelling the critical point values, as usual. 

It is not a difficult exercise to integrate equation ([19} 
and then transform the variable p in it to Cs, with the help 
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of equation (|21[). Once this has been done, the critical con- 
ditions, as given by equations p3|) . will have to be invoked, 
and all of these will deliver a relation for fixing the critical 
point coordinates in terms of the flow parameters, £ (which 
is actually Bernoulli's constant) and 7. This relation will 
look like 



7 + 1 



rcfl!)'(rc) + (f){rc) = £, 



(24) 



with the ranges of values of £ and 7 here, in what is es- 
sent ially a non-relativist ic approach, being accordingly cho- 
sen (|Das fc Sarkaij|200ll ). as opposed to the relativistic val- 
ues of £ and 7 adopted in Section [S] 

The choice of the pseudo-Newtonian potential, <j>{r), will 
obviously determine the number of roots of equation p4)) . 
Four such potentials have been considered here, and they 
have in general been labelled as = 0i(r), with {i = 
1,2,3,4}. In an explicit form, each of these potentials will 
be given as 



02 (r) 

03 (r) 
Mr) 



1 



2(r-l) 
1 



2r 



= -1 



In 



1-i 
r 



1/2 



('-;) 



(25) 



in all of which, the length of the radial coordinate, r, 
has been scaled in units of the Schwarzschild radius, de- 
fined as r-g = 2GM-BYi/(? ■ Every potential mentioned 
above has been introduced in accretion literature at vari- 
ous st ages to meet some specifi c physi cal requirement — 
1 by iPaczvnski fc Wiital lll980l). (bo. bvlNowak fc Wagonerl 
199ll ). and (^3 and 04 bv lArtemova et alj l lQQsh . respec- 
tively. With respect to spherically symmetric flows in par- 
ticular, a comparative overview of the phys ical properties of 
these potentials has been given bv lPas fc S arkar (2001) ■ 

Considering each of the potentials separately in equa- 
tion (|24|) . it will be seen that two distinct roots will be ob- 
tained on using both and (jia, while from (/)2 three roots 
will be delivered. The fourth potential, 04, will lead to a 
transcendental equation, and any root, therefore, can only 
be extracted by numerical methods. Using the bisection al- 
gorithm, it can be shown that only one physical root is pos- 
sible. 

With each such physically feasible root, a critical point 
can evidently be associated. The way to have any apprecia- 
tion of the behaviour of these critical points has already been 
outlined in Section [31 The first task would be to set up an 
autonomous dynamical system (in terms of a mathematical 
parameter, t), which will be 

d7(") 



= 2v^ let 



r0'(r)] 



dr / 

d7 = 



(26) 



Subject to the perturbation scheme. 



Csc + (Jcs and r — rc -I- (5r, equation (|26|) will lead to a set 
of coupled autonomous linear equations in the perturbed 
quantities hv^ and (5r, with having first been expressed 
in terms of iv^ and hr from the continuity condition, as 




Figure 3. Dependence of (always positive) on the parameters 
£. and 7 for the pseudo-Newtonian potential, 0i. The saddle-type 
behaviour is maximum for intermediate values of E and 7. 



„2 



1 

'2n 



(27) 



The coupled linear dynamical system will be 



with 



2c" 



2vtV Sr 



\ 2nJ n 



(28) 



V 



4c" 



From here it is an easy passage to deriving the eigenval- 
ues of the stability matrix associated with the critical points. 
With the use of solutions, Sv'^ ~ exp(f2T) and Sr ~ exp(57r), 
in equations (I28|l . these eigenvalues will be derived as 

i2 



c0'(rc)l^ 



(3 - 57) - (7 + 1) r< 



0"(rc) 



(29) 



from whose structure it can once again be claimed that the 
critical points can only be either saddle points or centre-type 
points. The dependence of Q,^ on £ and 7 has been separately 
shown in Figs. [3]and|4l under the choice of 0i and 04, re- 
spectively. For 01 there are two critical points, of which only 
one is the physically relevant saddle point, while for 04 there 
is only one critical point, which has to be found numerically. 
It is always a saddle point. The two potentials, 0i and 04, 
have been chosen because they give a closer Newtonian ap- 
proximation to fully general rela tivistic conditions, th an the 
other two potentials, 02 and 03 l|Das fc Sarkar|[200ll ). From 
both the plots in Figs. O and [IJ it is quite evident that con- 
trary to what it was for the fully general relativistic case, for 
pseudo-Schwarzschild fiows, strongly transonic features (in- 
dicated by high positive values of f2^) occur at much greater 
values of £ and 7. 



5 CONCLUDING REMARKS 

In trying to mathematically understand the nature of the 
three critical points in the general relativistic flow, it has 
been shown that there arises a situation, whereby, because 
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Figure 4. Dependence of (always positive) on the parameters 
£ and 7 for the pseudo-Newtonian potential, <j>4,. The saddle-type 
behaviour has a monotonic growth with increasingly higher values 
of £ and 7. 



of the fluctuating nature of the middle critical point (some- 
times a saddle point and sometimes a centre- type point), 
two contiguous critical points will be of the same kind. If 
these two particular points are both saddle-type then as an 
exercise of mathematical interest (if not of any direct phys- 
ical relevance) , it will be patently impossible to connect the 
two points by contirmous solutions. However, going back to 
equation (|8}, and subjecting it to a closer examination, a 
possible way of bypassing this difficulty could be found. It 
has been discussed already that the only critical conditions 
selected from equation (|8]) will be the ones on length scales 
greater than that of the event horizon. On the other hand, if 
some of the unphysical criteria {v^ = 1 at r = 1, or ?; = at 
r = 0, etc.) for criticality in the flow were to be taken into ac- 
count, then some new critical points (of mathematical inter- 
est only) could be found in the global phase portrait, even on 
unrealistic length scales. These possible critical points may 
then settle the difficulty which has arisen from the appar- 
ent existence of adjacent critical points of the same nature. 
Solutions could then be connected from one region to the 
other through these "hidden" critical points. 

Another problem with saddle points is that solutions 
passing through them are notoriously sensitive to the 
fine tuning of the outer boundary condition of the flow. 
This is a standing problem with stationary flows and in 
consequence of this, it has been shown for the classical 
Bondi problem that transonicity could be achieved only 
if the evolution of the flow wer e to be followed through 
time (|Rav fc Bhattacharieell2002h . This is a relatively easy 
proposition in the Newtonian domain. When a flow is stud- 
ied in the general relativistic regime, the time-dependent 
evolution will require much greater mathematical (both an- 
alytical and numerical) sophistication. Having made this 
point it should also be a fair expectation that transonic- 
ity would continue all the same to hold its primary position 
in spherically symmetric flows. 
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